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The rigid body motion of a two-body satellite under the oxtion of gravita- 
tional torques is analyzed. The satellite consists of two rigid bodies connected 
by a universal joint where damping is provided in the two journals. The 
motion of the satellite relative to the mass center thus has five degrees of 
freedom, two of which are provided with energy dissipation. It appears that 
the rigid body motion of such a composite satellite will automatically con- 
verge upon a motion in which a given axis of the satellite is earth-pointing. 

The equations of motion are derived directly from those of Eider. Neces- 
sary stability criteria are established. Numerical solutions for a practical 
scheme are presented. 

I. INTRODUCTION 

This paper deals with the analysis of the rotational motion of a satel- 
lite consisting of two rigid bodies connected by a hinge mechanism of 
universal joint type. The rotational motion of the satellite thus has five 
degrees of freedom; the two degrees of freedom that involve the relative 
motion between the two bodies arc provided with energy dissipation. 
It is found that any motion of the satellite with respect to the local 
vertical always involves relative motion between the two bodies. There- 
fore, the damping at the hinge joint dissipates not only the relative 
motion of the two bodies but also the motion of the satellite with respect 
to the local vertical. The satellite will then converge upon a stable mo- 
tion in which a specified axis of the satellite will remain close to the local 
vertical. 

The equations of motion are derived directly from those of Newton 
and Euler. This approach naturally suggests several additional de- 
pendent variables and results in numerically workable equations. This 
is not the case in the Lagrangian formulation. 

There are several practical problems involved in this scheme of pas- 
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sive gravitational orientation. One problem is to make the gravitational 
torque dominate over all other disturbing torques. A novel solution to 
this problem, which employs extensible rods, has been given by Kamm. 1 
Another problem is the development of the hinge dissipative mechanism. 
A viscous mechanism is described by Kamm, 1 whereas a hysteresis 
mechanism is suggested in this paper. These practical matters are not 
the substance of this paper; they are used as illustrations for the nu- 
merical treatment of a practical design. 

II. GENERAL EQUATIONS OF MOTION 

Consider a satellite which is constructed of two rigid bodies, with 
masses rai and m 2 , hinged at a point //. The centers of mass of the two 
bodies are denoted Si and S 2 , and the center of mass of the composite 
satellite is denoted S • Let the earth's center be and let Pi and P 2 
be arbitrary points of body 1 and 2. Also, denote OPi = Ri , OSi = 
Pi , OS = 9 , OS 2 = 92 , 0P 2 = R 2 , S 1 P 1 = ri , S 2 P 2 = r 2 , H& = £1 , 
HS 2 = £2 (see Fig. 1). (Note: £1 and £ 2 represent vectors, while A and 
(■1 which appear later represent their respective magnitudes. See Ap- 
pendix for list of symbols.) 

Let us introduce the following notations : 

"i > wii = angular velocity of body 1, 2, 

Th = reactive torque transmitted through the joint on body 1, 
F w = reactive force transmitted through the joint on body 1, 
Ti , T 2 = resultant torque on body 1, 2 exclusive of T w , 
Fi , F 2 = resultant force on body 1, 2 exclusive of F w , 
nij , m 2 = mass of body 1,2, 

■if) = m 1 ?n 2 /(wi -\- m->) = reduced mass of the system, 
m = mi + m 2 = total mass of the system, 
4»i , 4>2 = moment of inertia dyadic of body 1 , 2. 
Newton's and Euler's equations can now be written as 

Fi + F w = m/pi (la) 

F,. - F„ = m 2 Q 2 (lb) 

4>rwi + o)i X <IVo>i = Ti + T„ - £, X F„ (lc) 

4>,-oin + (on X 4» 2 -o>ii = T 2 - T// + £. X F„ (Id) 
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where the dots indicate time derivatives with respect to an inertial 
frame. Because of the constraint imposed by the hinge, the following 
relations are satisfied 

Pl = e + (£i - £ 2 ) — (2a) 

m 

or 

92 = 9 + (£, - £1) — . (2b) 

m 

Addition of (la) and (lb) yields the following vector equation which 
governs the motion of the mass center »S' : 

F, 4- F 2 = mi. (3) 

Using (la), (2a), and (3) we may solve for F w 

F w = ^F 8 -^F 1 + m(iJ 1 -£ 2 ). (4) 

in m 

Inserting (4) in (lc) and (Id) and using the fact that 

£i = g>i X £i and £ t = wi X £i + an X (g>i X £i) , 

etc., equations (lc) and (Id) become 

«I»/g>i + wi X */-u>i = T, 4- T„ 

+ £,X I— F, - - F 2 + m\t*n X (o>n X£->) 4- (in X £■•]> (5a) 
[ m in ) 

&2-am + unX*! ton = Tj — T,/ 

4- £ 2 X (- F 2 - - 2 F, + m[an X (a* X £i) 4- an X £.]) < 5b > 
where */ = 4», + fh{ttl — £,£.), i = 1, 2, and I is the unit dyadic. 

III. GRAVITATIONAL FORCE 

The earth's gravitational field is taken to be radially symmetric. The 
gravitational force, dG, , acting on an infinitesimal mass f/ra,- at P, is 

then 

,/G,= -^R, (6) 

where n = gR E ~ with y being the gravitational acceleration at the earth's 
surface and R E being the earth's radius. From Fig. 1 
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Fig. 1 — Vector displacement diagram of a two-body satellite. 
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where the last quantity represents terms of order t /pf and higher and 
I is the maximum linear dimension of the satellite. These higher-order 
terms are neglected in the analysis. Since Si is the center of mass of 
body % 

I r, dnii = 0. 

Hence Gi , the gravitational force on body 1, is 



pi" 



or, by (2a) 



Gi = |"_^ + ,_m (£a _ £i) . (I _ m J t + Q gyj (ga) 
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Similarly, 

G , = [-^ - •* (a - 4) ■ (I - 3W)][l + g)] (8b) 

where the symbol "*" denotes a unit vector. Using (7), the gravita- 
tional torque acting on body i about the center of mass is given by 



«- /r. X*,- J* X •,■*[! + og)] 



= 1,2. (9) 



Let F, = P/ + G, , T, = T/ + T Gi ,i = 1, 2. Substituting in (5a,b) 
with the gravitational torques in (9) and the gravitational forces in (8) 
with terms of 0(1/ p) and 0(f/p 2 ) neglected, the general equations of 
rotational motion of two hinged-connected rigid bodies become 

<!>/ • o>x + b)i X */ • o>i = -=■ p X */ • p 

P 

+ ^(£i X £2 - 3£, X pp-£ 2 ) +T/ + TH 

P 1 (10a) 

-5 il xF,' + 5«F 1 ' 

+ m(fa>n-£2<£i X on — w n 2 £i X £2 + J2i'£a«n — Wrwn), 
*2' • (i>n + "n X W • wii = — p X 4»2 • p 

+ ^ (£ 2 X £1 - 3£ 2 X pp-£i) + T,' - T„ 

P 3 (10b) 

+ m(wi-£i£2 X o>i - w^JBi X £1 + £ 2 -£iu>i - £i£ 2 -wi). 

Note that TV and T 2 ' are the resultant torques imposed on body 1 and 2 
by some external sources other than gravity. They do not include 
torques arising from the reaction of one body upon the other. Simi- 
larly, F/ and F 2 ' are the resultant forces on body 1 and 2 due to ex- 
ternal sources other than gravity. They do not include the reaction of 
one body upon the other. Thus both these torques and forces are worked 
out as though the bodies were not connected. Various environmental 
disturbances, like solar radiation pressure or interaction of a magnetic 
moment in the satellite with the geomagnetic field, may be taken into 
account by assigning appropriate values to T/, T.', F/, and F>'. This 
subject is not treated here to conserve space. 
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From (8a) and (8b) it is seen that 

Gl+ G 2 = -«. (lla) 

p 3 

If the gravitational forces are the only ones in Fi and F 2 , then (3) 
becomes 



.. _ MP 



(Hb) 



The solution of this vector equation is an elliptical orbit of S , inde- 
pendent of the rotations of the satellite, because of the fact that terms 
of 0(l 2 /p 2 ) have now been neglected. As the earth's gravitational field 
is assumed to be radially symmetric, the orbital plane is fixed in the 
inertial space. 

IV. COORDINATE SYSTEMS 

Four reference frames are used to describe the motions of the satellite. 

The first frame has its origin at the geocenter with the Z-axis 
through the perigee of the orbit and with the F-axis in the direction of 
the orbital angular momentum. The X-axis is chosen to form a right- 
handed set of axes (see Fig. 2). This coordinate system is taken to be 
inertial. 

The second is an earth-pointing frame. It has its origin at the satel- 
lite's center of mass, S , with the z-axis along OS making an angle \f/ 
with the Z-axis. The y-axis is parallel to the F-axis, and the re-axis is 
chosen to form a right-handed system. The relationship between the unit 
vectors of the coordinate systems O-XYZ and So-xyz is 

/x\ (c+ o -sA(x\ 

[y - 1 F (12) 

w w o WW 

where S and C are abbreviations of sine and cosine. 

The third frame has its origin at Si with axes Si-XiyiZi along the 
principal axes of inertia of body 1. Euler parameters are employed to 
describe the motion of Si-XiyiZi relative to Sg-xyz. The transformation 
is given by 



(13a) 
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Fig. 2 — Coordinates of the rotating and nonrotating frames. 



where 



t* _2 _ 



(an) = ( 2(fc, - fx) 
<2(& + vx) 



f 2 + X 2 



2(fe + fx) 

-$ 2 + »? 2 - f 2 + X 2 
2(-£x + tf) 



2(tf - ux) 
2(fx + tf) 

_£2 _ ,2 _J_ f2 + x 



.) 



(13b) 



i, j = 1, 2, 3 representing rows and columns respectively, and 

s 2 + rf + r 2 + x 2 = 1. (13c) 

The fourth frame has an origin at S 2 with axes Sz-wjiZo along the 
principal axes of inertia of body 2. If a universal joint is used, the relative 
rotation of the second body can be completely specified with only two 
angles, namely a, the rotation of the journal in body 1, and /3, the rota- 
tion of the journal of body 2. When these two journals are directed re- 
spectively along Xi and t/ 2 , then the transformation from Si-XifjiZi to 
SrXiV&i is given by 



= (bn) 



(14a) 



where 



(M = 



Sa Sfi 
Ca 

-SaCP 




(14b) 



The constraint equation Xi-y 2 = is automatically satisfied by the in- 
troduction of the two coordinate parameters a and /3. The angular 
velocities of the two bodies are 



*■>! = H/ + ^'i + ^2/1 + A3Z1 



(loa) 



2246 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1963 

where 

Xi = 2( X £ + N - vt ~ &) (15b) 

A 2 = 2(-f| + X »» + tf ~ »?*) (15c) 

X 3 = 2(^ - ft + xf - rx) (15d) 



and 



ton = "i + ah + #/o . (15e) 



V. SPECIALIZED EQUATIONS OF MOTION 

Let us specialize our satellite so that £2 = and £1 = — (fa . We 
assume gravity to be the only external force, i.e., T/, T 2 ', F/ and F 2 ' in 
(10) are taken to be zero. Then equations (10) are equivalent to those 
derived from two bodies connected at their centers of mass except that 
the inertia dyadic <I»i is replaced by * x ' denned in (5) (*&•/ = *2 as 
(•> = 0). The two bodies are connected by a universal joint, which is 
characterized by an interposed weightless body, having two perpen- 
dicular journals as previously described. The torque T H , transmitted 
through the universal joint, consists of the constraint torque T* , the 
elastic restoring torque T r , and the dissipative torque T^ . The com- 
ponents of the latter two along the journals X\ and y 2 are specified by 
subscripts 1 and 2 respectively. Hence T H can be written as 



T„ = TA X fa + (T rl + Tufa + (T r2 + T d2 )fa . (1(5) 



Let 



/i = *i'-& (17a) 

h = #i'-fr (17b) 

h = &x'-zi (17c) 

7.i = 2 -x 2 (17d) 

h = * 2 -y 2 (17e) 

h = &2-Z-2 (17f) 

Wi(i = 1, 2, 3) = components of toi along :i\ ,y\,Z\ (17g) 

a-,(t = 4, 5, 6) = components of o>n along f 2 , ^ 2 , z* . (17h) 
From the orbit equation (lib), the following relations can be derived 

* = n ° o V . (1 + «W 2 (18) 
(1 - e z ) s 
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3m _ 3fl 2 

P 3 (1 - <?Y 



G = ^ = TT ™— ■ (1 + ^)' (19) 



where 



e = eccentricity of the orbit 

tt = 2ir divided by the orbital period. 

Euler's equations of motion (10), simplified for the specialized satellite, 
are written out as 

/.ci, = (/, _ / :( ) (tojjw,, - Gwiz) + T Ti + Tat (20a) 

/,u 2 = (/., - /,) ( W3W , - r/// ;t / h ) + (T„ + 7V-X '« - T e Sa (20b) 

fab = (/, - h) (wm - Qthih) + (Trz + Z\«)fla + T c Ca (20c) 

/ 4 ci 4 = (h - h) (wewe - On**) - (r„ + I\a)C/3 + T C SP (20d) 

/ 5 ^5 = (/c - A) Uuu - Gn#u) - Ta - T d2 (20e) 

/ 6C i 6 = (/ 4 _ / 5 ) ( um - (Jri^) - (T rl + T A )8fi - T e Cp (20f) 
where 

«,• = o,-3 (see 13b) i = 1, 2, 3 (20g) 

n,^ = E 3 Mti (^e 14b) t = 1, 2, 3. (20h) 

fc-I 

Because of the constraint £•# = 0, a relation must exist among the 
six oj.'s. Such a relation, i.e., (g>i - u n )(£i X 2/2) = 0, can be ob- 
tained from (15e). This yields the following relationship: 

oi,Sa - u 3 Ca - mS0 + w«C/3 = 0. (21) 

If (21) is differentiated and equations (20) are substituted, the un- 
known T, is found to be 



._ (fa e* , sv cy\ 



Y K/ 3 - /i)(«i«a - GW*) + (^n + T d2 ) Ca] 
_ ^ [(/, - 7,)(«i«i - Gnin 2 ) + (ZV2 + r«) &*] 

is 

_ #* [(/ 6 _ / 6 )(« B «« - (The*) - (Tri + T dl ) C{3] 
1a 
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+ £? [(/ 4 - / 5 )Uu 5 - Gn 4 n b ) - (7' rI + T dl ) S0\ 

+ a{^C<x + o> 3 &*) - ^(co 4 C|S + o. 6 S/3) > . (22) 

Equations (20) could be considered as a system of six second-order 
equations in six unknowns £, -q, f, x, «, /3, while ^ is determined from (18) 
and the w»-'s from (15). For computation purposes it is convenient to 
also leave the o>,'s as dependent variables. Equations (15) give 

£ = KxXi - fX, + >?X») (23a) 

■h = i(fXi + XA* - frs) (2.3b) 

£ = *(— ifXi + £X 2 + X a 3 ) (23c) 

X - *(-£X, - rjXo - fX,) (23d) 

a = -an + W4 Cj8 + o> fi tf/3 (23e) 

# = — cooCa — a>3»S'a + Wft (23f) 

X,- ■ «< - aal, i = 1, 2, 3. (23g) 

There are now 12 first-order equations, (20 a-f) and (23a-f), in the un- 
knowns £, r/, f , x, «, 0, a-i ,0)2,0)3,0)4,0)5, and o) C . If Euler angles had 
been used instead of Enler parameters, there would be certain positions 
of the body for which the derivatives of the angles have a singularity. 
However, no singularities occur when Euler parameters are used, as can 
be seen from (23). It should also be noticed from ( 13b) that the matrix- 
fixing the body position is not changed if the coordinates (£, rj, f, x) are 
replaced by (-£, -77, -f, - x ). 

VI. DISSIPATIVE AND ELASTIC TORQUES IN THE UNIVERSAL JOINT 

To completely define the problem it is necessary to specify the elastic 
and dissipative torques T r and Td . 

6.1 Damping Torques 

Two types of damping torques are considered here. The first is viscous 
damping of the linear velocity type; the torque on body 1 has two com- 
ponents 

T dl = CMi (24a) 

T d2 = C 2 /3& (2 IM 

where C\ and C 2 are viscous damping coefficients. 
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The second is of magnetic hysteresis type. The damping is furnished by 
hysteresis losses produced by the relative motion of a permanent mag- 
net and a permeable material. The torque in the ^-direction might be 
approximately expressed by the following process. If a > 0, the torque 
would be represented in region I in Fig. 3 by 



r - t * J- T a ~ a 
1 «ji = I <fi -r * di z — 



(25) 



as long as | T dx \ < f d i where a, f d \ are constants and a*, T d * are the 
values of a, T dx when a last changed sign. After | T di | reaches 7' dl then 
T d i remains at f dl as long as a does not change sign. This is represented 
as region II in Fig. 3. If a changes sign, then (25) applies and the proc- 
ess is repeated. This is represented by region III of Fig. 3. T d2 is defined 
by replacing a by /3 and subscript 1 by 2 in (25). According to this 
idealized hysteresis, no energy is dissipated in region III. In an actual 
device, energy would also be dissipated in this region because of minor 
hysteresis loops. The chief advantage of magnetic hysteresis damping 
is that it is amplitude dependent instead of velocity dependent, since 
the librational frequency, which is of the order of the orbital frequency, 
is too low to make the velocity damping effective. Other merits of the 
magnetic hysteresis damper will be stated in the descriptions of a prac- 
tical design for a numerical computation. 



T d , 




*-a 



Fig- .'5 — Magnetic hysteresis damping torque produced by a magnetic device 
on the aii- journal. 
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6.2 Elastic Torques 

It is assumed that each journal is furnished with a linearly elastic re- 
storing torque produced, for example, by the torsion of a wire. The torque 
acting on body 1 is given by 

T ;1 = hah (26a) 

T, 2 = fctffr . (26b) 

where &i and k 2 are spring constants. 

VII. MISCELLANEOUS TORQUES 

Many other torques such as those due to interaction of the satellite's 
magnetic moments with the geomagnetic field, solar radiation, self- 
gravitation between two bodies, and plasma effects will act as forcing 
terms in the equations of motion. By proper design, these torques can 
be made small compared to the gravitational torque. However, since the 
gravitational torque varies inversely as the cube of the geocentric dis- 
tance, it may not necessarily dominate in the orientation of satellites in 
very high orbits. Also, in very low orbits, aerodynamic drag may be big 
enough to upset the orientation. If long rods are used with weights on 
the ends, the gravitational torque can be made to dominate for a certain 
range in altitude. 

VIII. EQUILIBRIUM AND STABILITY 

Let us consider only the equilibrium position 

(^,7?,f,X,«,/3) = (0,0,0,1,0,0) 

ill which the Xj , ij\ , z x axes are lined up with the x,y,z axes. For viscous 
damping, the stability criteria for the position (0,0,0,1,0,0) can be found 
by linearizing the equations of motion about this position. The same sta- 
bility criteria are obtained for equilibrium positions found by rotations 
of 180° around the x, y, and z axes, i.e., (1,0,0,0,0,0), (0,1,0,0,0,0), 
(0,0,1,0,0,0). For hysteresis damping, there will be an infinite number of 
stable equilibrium positions. All of these can, however, be made suffi- 
ciently close together to either one of the above four equilibrium posi- 
tions, thus maintaining an axis in the satellite nearly in line with the 
local vertical. 

From the definition of Euler parameters, the infinitesimal angles of 
rotation about the X\ , //i , and Z\ axes are £i = 2£, 771 = 2?;, ft = 2f, de- 
fined as the roll, pitch and yaw angles. If £•_> , rj- 2 , ft are the infinitesimal 
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angles that the principal axes of body 2 make with respect to the rotat- 
ing coordinate system S Q -.ryz, and a and /3 are small, then 

£2 = a + & (27a) 

772 = + Th (27b) 

f 2 = fj ■ (27c) 

In the linearization process, we take the eccentricity of the orbit, e' 
to be small in order to insure the realization of the infinitesimal angles' 
This is necessary in view of the well known result of the satellite pitch 
motion that the angular excursion produced by the eccentricity is of the 
same order of magnitude as the eccentricity itself. From (18) ty becomes, 
with zero phase angle, 

yf, = m + 2cSm + 0(e 2 ). (28) 

To linearize the general equations of motion given by (10), let us assume 
viscous damping as expressed in (24) and linear restoring torques as 
given in (26). The perturbing torques and forces, TV and F/ (i = 1, 2), 
are neglected. Also, let £1 = -/A and£ 2 = 4h . Then, equations (10) 
are linearized to the following: 

m + lam + &'(* ~ fr) + ** " *i'l" = 2<tf < 1 + Ll)Sm < 29a > 
fa + im + <V(* - *) + A* - **'*» - 2etf ( 1 + k)** ( 29b ) 
£1 + 2Vifc + Ci"(ki ~ kt) + ( M + wife - £& = (29c) 

h + tfJh + &*(& - &) + flrftf + "* " &6 = (29d) 

f + (1 - A - JV)"T - n/ifc - D/* = ° ( 29e ) 

where 
Jj X = ifi (J.,/ L, , In = 'Tidfo/h 

Ci = C-i/h, Co' = Ct/U 

W = h/h , hi = h/h 

d x = 3sr(/i - D/h + 30% + hi 
<h = 3<f(/ 4 - h)/h + 30% + KV 
AT! = /»r/o//, , A~ 2 = WV2//4 

c/' - c,//i , f 2 " = cy/ 4 
r/, = (Ji + h - h)/h 
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q-2 = (h + h - h)/h 

w, = 4fi 2 (l - g, + 3A/,/4) + AV/i 

Uj - 4S2 2 (1 - q 2 + 3A/ 2 /4) + fci/*« 

£, = fci/7, - QTNi , £ 2 = k/J 4 - S2 2 JV 2 

/i = (/i + h - h)/(h + h), /, = (/ 4 + J. - /,)/(/, + / 6 ). 

It should be noticed that the pitch equations (29a,b) do not depend 
on & , £ 2 , f and are decoupled from the roll and yaw equations (29c,d,e). 
The eccentricity enters as the amplitude of a forcing term in pitch but 
not in roll and yaw. The transient part of the pitch libration can be solved 
from (29a,b), excluding the forcing terms, by substituting with 

77, = Bifi", i = 1,2. 

The resulting characteristic equation in s is then 

(1 - UL^s" + (C7 + CV + Cx'U + C'LJs* 

+ (rfi + d 2 + h'U + Th'Li)i 

+ (d x C/ + rf 2 C/ - CW - C 2 %')s + (did* - h'h') = 0. (30) 

The pitch motion is damped about (0,0,0,1,0,0) if and only if the Routh- 
Hurwitz conditions are satisfied. This insures that the real parts of the 
roots of (30), representing the damping constants for the two principal 
modes, are negative. These give 



where 



h > h 


(31a) 


hi > — 312 - - (/ 4 — h + mhh) 

1 x 1 z 


(31b) 


h — h + wiA^ ^ 7 4 — /e + *MiA 


(31c) 


/, = J x + I 4 + 2A^ S 




J, = h + 7 6 + 2mf 1 ( 2 




L = h + h. 





Ix , Iy , Jz represent the moments of inertia of the composite body about 
*S'o . Condition (31a) is the same as that of a single rigid body. Condition 
(31b) states that /»•> must be larger than a certain critical value if one 
body is unstable (e.g., h — I 6 + fhtiU < 0). This value is zero if both 
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bodies are stable by themselves. It can be shown that there cannot exist 
a cocked equilibrium position in pitch if the parameters are such as to 
make the position (0,0,0,1,0,0) stable. Condition (31c) implies that 
there exists an undamped motion if the equality sign holds. This rigid 
body motion has a frequency w, , given by 



= 39: 



(/i - h + mdh) 



(32) 



I 2 + mt^ 

The roll and yaw equations (29c,d,e) are all coupled. This justifies 
the use of a damper only for roll. Up to first-order terms, there are no 
forcing terms due to eccentricity. The characteristic equation is 

/; 6 s 6 + 6 5 s 5 + 6 4 s 4 + &*' + b 2 s 2 + b v s + b = (33) 

where 

b Q = tt 2 (l -/, - f. 2 ) ( Ul n 2 - kik 2 ) 

6, = i2 2 (l - h - ft) (C2V + C\"u 2 - hCY - hC,") 

bi = ft 2 (l - /i - ft) (ui + w 2 + Nik 2 + NA) + wi«2 - hh 

+ tf(fiqik\ + hqjii + /tfi^ + /ifttO 
h 3 = « 2 (1 - f v - f 2 ) (&* 4- r 2 " + C 2 "AT, + C,"AT 2 ) 

4- &"(«» + tf/172 + tf/«g» - &) 

+ C 2 "(ui + flVtfi + tffifc - &) 
64 = fl 2 (l - /, - /«) (1 - N1N2) + Mi + u 2 + N& + iV 2 fci 

+ rf(/ift + ./Vy 2 - /i^Vi - faiNi) 
b b = CV'(1 + N 2 ) + C/(l + AM 
6 fl = 1 - NiN 2 . 
The Routh-Hurwitz stability criteria are 



bo > 0, bx > 0, 
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b fi 
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> 0, 6 6 > 0. 



b 6 b 5 

If these are satisfied, there will be three more modes of damped libra- 
tions. Due to the coupling between the roll and yaw librations, the yaw 
libration can be damped out by the roll damping, as can be observed 
from (29c,d,e), although no yaw damping mechanism is provided in the 
present scheme. Hence, all modes can be damped out and the satellite 
will oscillate with some steady-state amplitude about an equilibrium 




Fig. 4 — Gravitationally oriented two-body satellite with extensible rods. 
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position. Some of these conditions are too complicated to give any physi- 
cal insight. However, some are quite simple and are given below. 

Since the parts of b , b 2 and 6 4 which involve h are (fa/C^bi , 
(/ci/Ci)& 3 , (fci/Ci)&5 respectively, multiplying the odd columns of the 
Hurwitz determinants by h/Ci and adding to adjacent columns will 
eliminate the fcj terms. Hence the only condition on fci is b > 0, i.e., 



fci > 12' 



[4(/ 2 - h + f md(2)(h ~ h + f mA/s) - | mtiti] 
L - h 



(34) 



As fci and /c 2 approach infinity, the satellite becomes one rigid body. Since 
the stability conditions are not changed by an increase of fci (and k 2 ), 
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Fig. 5 — Angular variation between the 2i-axis of the satellite and the local 
vertical for a hysteresis damper, cos 0. 
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it appears that the single rigid body criteria for roll and yaw stability 
are necessary. These are 



where 



P X P Z > 


(35a) 


1 + SP X + PxP Z > WPrPz 


(35b) 


'u ** p *v 1* 





/ 



L 



and J i , I y , and 1 z are given in (31). Condition (35a) can be verified 
from the inequality h\ > 0. Other necessary conditions in the case of 
4 = are found from the third-order Hnrwitz determinant to be 
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Fig. 6 — Relative angle about the rci-journal for a hysteresis damper, 
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IX. BISTABILITY 
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(35c) 

(35d) 



The satellite is in a stable equilibrium position if the z^axis is in line 
with the local vertical (i.e., the z-axis) pointing in either direction. If a 
directional device such as an antenna or a camera is used along the 
negative Zi-axis, it may point at or away from the earth. The equi- 
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Fig. 7 — Relative angle about the j/ 2 - journal for a hysteresis damper, /3. 
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librium positions (0,0,0,1,0,0) and (0,0,1,0,0,0) correspond to the device 
pointing toward the earth, whereas (1,0,0,0,0,0) and (0,1,0,0,0,0) cor- 
respond to the device pointing away from the earth. In the latter case 
an inertia wheel in the satellite can be activated with a predetermined 
number of turns, and the satellite can be rotated 180 degrees so that the 
device will be earth-pointing. The equations governing this turning are 
given by (10), where the applied torque on body 1 is approximately 



T/ = -^-[JMSxt + SSy,)} 
at 



(36) 



where J m is the angular momentum of the inertia wheel and 8 is the angle 
between the £i-axis and the axis of the inertia wheel. Another scheme 




Fig. 8 — Component of angular velocity of the satellite along the yi-axis for 
a hysteresis damper, w 2 /12. 
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would be to use two devices, one on each side of the satellite, directed 
along the positive and the negative direction of the Zi-axis respectively. 
Only the one that is earth-pointing would be activated. 

X. NUMERICAL RESULTS OF A PRACTICAL SCHEME 

A practical scheme, as shown in Fig. 4, is suggested here for a com- 
munications satellite. The particular construction, employing extensible 
rods and tip masses, is to effect large moments of inertia so that the 
gravitational torque will dominate over all disturbing torques. Body 1 
of the satellite, which consists of the satellite's main structure (with 
directional antennas) and a mast rod, is to be earth-pointing. Body 2, 
being an auxiliary body for attitude-control purpose only, is constructed 
of two rods and is in an unstable position with respect to the local vertical. 
These rods arc extended, upon ejection from the launching vehicle's 
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Fig. fl — Angular variation between the 2i-axis of the satellite and the local 
vertical for a viscous clamper, cos 6. 
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final stage, by unrolling from sheet metal drums. The universal joint 
employs torsion wires to produce elastic restoring torques and provides 
hysteresis damping by relative displacement between magnets and a 
permeable material. (See Fig. 5 of companion paper. 4 ) 

The advantages of magnetic hysteresis damping are that it is ampli- 
tude dependent, insensitive to temperature variation, involves no sliding 
parts and requires little weight. Coulomb friction damping, while also 
amplitude dependent, is less desirable because of possible cold welding of 
sliding parts in the high vacuum of space. Velocity-dependent damping 
by employing viscous fluids is believed to provide lower damping for a 
given weight, and the viscous fluids involve questions of temperature 
sensitivity. 

All the parameters are chosen based on the adjusted moment of inertia, 
1 1 , of body 1 subject to stability criteria and other necessary considera- 
tions. The stability criteria (31b) and (34) specifying the critical values 
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Fig. 10 — Relative angle about the zi-journal for a viscous damper, a. 
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Fig. 11 — Relative angle about the j/s-journal for a viscous damper, 0. 






of lc-i and k\ , respectively, which are derived from viscous damping, are 
found to apply approximately also in the case of hysteresis damping. 
These parameters are: 7,//i = 1.00, 0.003, 0.159, 0.381, 0.540 
(i = 2, ■■■ , fi); A-,//iS2 2 = 1.131, 2.238 (* = 1,2); for a hysteresis 
damper: T di /h& = 0.159, 0.210 (i = 1,2), a = /3 = 2°; for a viscous 
damper: Ci/IiQ = 0.870, 1.281 (i = 1,2). With the alcove value of 
the viscous constant C-> , the amplitude of the lower mode of pitch 
libration can be reduced according to (30) by a factor of e in 0.22 orbit, 
which is close to the optimum. The optimum in the case of pitch mo- 
tion was found by Zajac 5 to be 0.137 orbit. Equations of motion (20) f 
depend only on the above dimensionless parameters and are independent 
of I\ and £2 as long as / is measured in fractions of an orbital period. Some 
initial conditions which might simulate a micrometeoroid impact or the 



f Equations (20) were programmed on an IBM 7090 by Mrs. W. L. Mammel. 
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motion after the erection of the rods are at £ = : £ = 77 = f = a = /3 = 0, 
X = 1, «i ™ fi, ai2 = 5fi, 0)5 = ft, o> 3 = o>4 = co6 = 0. Figs. 5-8 represent 
the computer solution of equations (20) using a magnetic hysteresis 
damper. In Fig. 5, is the angle between the zi-axis and the local vertical. 
The satellite stops tumbling after four orbits and settles to within 10° of 
the local vertical after six orbits. The satellite librates about a cocked 
equilibrium position indefinitely due to the forcing torque of orbital ec- 
centricity (e = 0.01). The pitch angular speed of body 1, o> 2 , approaches 
one revolution per orbit, which is the proper speed for an earth-pointing 
satellite. Figs. 9-12 show similar results of a viscous damper. In this 
case the satellite ended up in an inverted position. 

Effects of the environmental disturbing torques, such as those due to 
solar radiation and the interaction of the magnetic moment in the satel- 
lite with the geomagnetic field, have been investigated, although the 
results are not included here. Cases with various other initial conditions 
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Fig. 12 — Component of angular velocity of the satellite along the f/i-axis for 
a viscous damper, w_>/£2. 



TWO-BODY SATELLITE 22G3 

have also been computed. All these results indicate that gravitational 
orientation of a two-body satellite is feasible. 

APPENDIX 

Nomenclature 

A.l Latin Symbols 

an = direction cosines of Si-XjyiZi frame with respect to S - 

xyz frame (i,j = 1,2,3) 
bn = direction cosines of Si-sty^i frame with respect to S\- 

x&iZi frame (i,j = 1,2,3) 
bi = coefficients of characteristic equation of £1 , £2 , f (i = 

0, 1, •••,0) 
Bi = complex constant of ??,• (* = 1,2) 
C = cosine operator 

C'i = viscous damping constants of a,P (i = 1,2) 
Ct, C" = adjusted damping constants of a,/3 denned in equations 
(29) (i = 1,2) 
di = coefficients denned in equations (29) (i = 1,2) 
ft = moment of inertia coefficients defined in equations (29) 

(i = 1,2) 
F w = force on body 1 due to reaction of hinge 
F,- = resultant force on body i exclusive of F* (i = 1,2) 
F,' = resultant force on body i exclusive of gravity and F fl 

(* = 1,2) 
g = acceleration of gravity on the earth's surface 
G = quantity denned in equation (19) 
G,- = gravitational force on body i (i = 1,2) 
H = hinge point 
I = unit dyadic 

/, = adjusted moments of inertia (i = 1, • • ■ ,0) 
J x , I y , I z = moments of inertia of composite body about the common 
center of mass 
./,„ = angular momentum of inertia wheel 
ki = spring constants producing torques in X\ , y-i directions 

(t = 1,2) 
h/,ki = adjusted spring constants denned in equations (29) 

(* - 1,2) 
I = maximum linear dimension of the satellite 
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£, = position vector of center of mass of body i from hinge 

(* - 1,2) 
(i = magnitude of £ t - (i = 1,2) 

Li = coefficients denned in equations (29) (i = 1,2) 
m,- = mass of body i (i = 1,2) 
m = total mass of satellite 
m = reduced mass 
iii = direction cosines of z-axis on Si-XiffiZi and SrXiy&i frames 

(i = 1, • • • , 6.) 
AT,- = coefficients denned in equations (29) (i = 1,2) 
= center of the earth 
Pi = arbitrary point in body i (i = 1,2) 
Pi , P v » P* ~ ratio of moments of inertia in equations (35) 

(/,• = coefficients denned in equations (29) (i = 1,2) 

Ti = position vector of Pi from center of mass of body i {i = 

1,2) t 
R, = position vector of P.- from (i = 1,2) 
R K = mean radius of the earth 
*S = sine operator 

s = variable in characteristic equations 
So = center of mass of satellite 
Si = center of mass of body i (i = 1,2) 

/ = time variable 
T„ = reaction torque transmitted through the joint on body 1 
T; = resultant torque on body i exclusive of T„ (i = 1,2) 
T/ = resultant torque on body i exclusive of T w and gravita- 
tional torque (i = 1,2) 
T,n = gravitational torque on body i (i = 1,2) 
T ( . = constraint torque of joint on body 1 
Ta = dissipative torque of joint on body 1 
Tj, = magnitude of saturated hysteresis torque of magnet i 

a = 1,2) 

Tdi* = value of Td-. when a (i = 1 ) and (i = 2) last changed 
sign 
T r = elastic restoring torque of joint on body 1 
u,- = coefficients defined in equations (29) (i = 1,2) 
X,Y,Z = fixed frame coordinates 
x,y,z = rotating frame coordinates 
Xj , i/i , Z\ = body 1 coordinates 
Xt , yi , z-2 = body 2 coordinates. 
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A.2 Greek Symbols 

a = relative angle of rotation of body 2 about £i-axis 

5 = constant of magnet 1 
a* = values of a when a last changed sign 

= relative angle of rotation of body 2 about y 2 -axis 

/3 = constant of magnet 2 
0* = values of when last changed sign 

5 = angle between avaxis and the inertia wheel axis 
e = eccentricity of the orbit 

f = Euler parameter 

f, = infinitesimal angle about 2,-axis (t = 1,2) 

■q = Euler parameter 

m = infinitesimal angle about i/.-axis (i = 1,2) 

6 = angle between z r axis and the local vertical or z-axis 

X j = components of the relative angular velocity of body 1 with 

respect to rotating frame (i = 1,2,3) 
M = a gravitational constant of the earth 
£ = Euler parameter 

£< = infinitesimal angle about a\-axis (i = 1,2) 
p = position vector of S from 
p, = position vector of Si from (i = 1,2) 
«I», = moment of inertia dyadic of body i (i = 1,2) 
<!>/ = quasi moment of inertia dyadic of body i (i = 1,2) 
X = Euler parameter 
\f/ = true anomaly of ellipse 
12 = mean orbital angular speed of satellite 
»i , «n ™ angular velocity of body 1 ,2 
coi , wi , oj :i = components of a>i along #i , y% , z\ axes 
a-4 , oifi , coo = components of un along x* , ],-> , z-> axes 

u), = natural frequency of an undamped roll libration. 



A.3 Notes 

" = unit vector 



(-D 



• = time derivative in an inertia! frame ( = 

boldface characters indicate tensors and vectors (it is assumed that 
dropping the boldface means the magnitude of the vector; i.e., p = 

I el). 
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